Meteoroid bulk density and Ceplecha types 
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Determination of asteroid bulk density is an important aspect of NEO 
characterization, yet difficult to measure. As a fraction of meteoroids originate from | 
asteroids (including some NEOs), a study of meteoroid bulk densities can 
potentially provide useful insights into the densities of NEOs and PHOs in lieu of 
mutual perturbations, satellite, or expensive spacecraft missions. 


NASAss Meteoroid Environment Office characterizes the meteoroid environment for | 
the purpose of spacecraft risk and operations. To accurately determine the risk, a 
distribution of meteoroid bulk densities are needed. This is not trivial to determine. 

If the particle survives to the ground the bulk density can be directly measured, 

_ however only the most dense particles land on the Earth. The next best approach is 
to model the meteor’s ablation, which is not straightforward. Clear deceleration is 

| necessary to do this and there are discrepancies in results between models. 


densities, then calibrate the proxy with known densities from meteorite falls, 
_ ablation modelling, and other sources. An obvious proxy choice is the Ceplecha 
type, Kz, thought to indicate the strength of a meteoroid. 


One approach to a distribution of bulk density is to use a measured proxy for the 


K, is frequented cited as a good proxy for meteoroid densities, but we find it is 
poorly correlated with density. However, a distinct split by dynamical type was seen 
with Jovian Tisserand parameter, T,, with meteoroids from Halley Type comets (T, 
< 2 ) exhibiting much lower densities than those originating from Jupiter Family 

\ comets and asteroids (T, > 2). 
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In 1958, Ceplecha introduced a parameter, K,, that he saw as a measure of the 
strength of a meteor and was linked to meteoroid densities (Ceplecha, 1958). 


It was noted that the beginning heights, H,, 
of meteoroids were correlated to their 
strength based on measured H, of more 
porous meteoroids (Draconids) or more 
dense meteoroids (Geminids). H, is most 
affected by velocity, but separate bands can 
clearly be seen when plotting H, against 
velocity. Ceplecha thought these bands were , 7 
correlated with meteoroid strength. esis 
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Figure 2: Heat map of beginning height vs. velocity for 
sporadic meteors seen in NASA's Wide-field camera 
network. 


What is K,? Based on a expression from Levin (1956) for surface temperature of a 


meteoroid per height: T(0,h) is surface temperature at a height 


A is heat conductivity 
6 is meteoroid density 
1 (OM imu elem-Leeoysstenteler-lulesmeelualeleels 
COS TAI VA R) c is the specific heat of the meteoroid material 
b is the air density gradient 
Zp is the zenith distance of the radiant 
V.,. is the approach velocity 
o is the air density 
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If the surface temperature and air density are set to values at the meteor 
beginning height. Ceplecha put all the physical constants on one side and 
observable quantities on the other, and set the physical constants are set to K,: 


Kp = log(pz) + 2.5 log(v.,) — 0.5log(cos(Zp )) 


Where K, is a function of material constants and surface temperature, and changes 
in K, are strongly tied into the composition of the meteor. 


Ceplecha (1958) found several K, groups that show up when isolating the meteors 
and K, values by velocity and orbital type, as seen in Table 1 and Figure 3. 


Table 1: Groups of meteors according to Ceplecha (1988). 


Ceplecha Group K, Range Orbital Criteria 


Asteroidal 8.00 < K, 
Group A 7.30 < Kz < 8.00 
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In order to use Ky as a proxy for density distribution, the link between K, and 
density must be established. The earliest studies on meteoroid bulk density were 
from the 1950’s. Since then, many studies have been performed. 


=O. 
Group B 7.10 < K, < 7.30 q < 0.30 AU 


Group Cl 
Group C2 
Group C3 
Group D 


6.60 < K, < 7.10 a <5 AU,i = 35° 

6.60 <= K, <= 7.10 a2=5AU 

6.60 <= K, <= 7.10 ad ——o AU, 1 > 35° 
Kz < 6.60 


Figure 3: Distribution of K, as seen in NASA 
Wide-field meteor cameras. | 


These groups show up clearly on a H, vs. velocity graph, further confirming 
Ceplecha’s suspicion they are tied to the meteoroid strength. 
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In examining the various sources that have linked density to K,, the relationship 
was not as clear as one would hope and primarily based on studies that have 

_ significant biases or models that have since been disproven, such as single-body 
ablation. Other discrepancies likely stem from a paucity of data. Several of 
Ceplecha's papers relied on meteorite density data from two or three meteorite 
falls, or just one fireball (Ceplecha, 1977, 1988; Ceplecha et al., 1998). Table 2 
shows seven of the best studies. Even within these studies, much is unclear and 
unknown. Note the discrepancies, particularly in Groups A through D. 


Table 2: Average meteoroid density corresponding to different K, groups. 


Density by Ceplecha Group (g cm?) 
Study DYeks inva) Com eveye! 


Ceplecha Rough estimates using past 


(1966, 1967) work 
Single-body ablation cto) 0.93 0.34 


Revelle (1983) 


Ceplecha Three meteorite falls, other 1.4-2.7 0.65-1.7 0.55-0.91 0.18-0.38 


(1993) theoretical assumptions. 
Ceplecha et al. 3 meteorite falls, known meteor 1 0.75 O27 
(1998) shower parent bodies 


Revelle (2001) Uniform bulk density model + 
inhomogeneous porous model 


0.75-0.93 0.27-0.34 


Bellot Rubio et 
al. (2002) 


Kikwaya et al. Ablation modelling with ORs} 
(2011) thermal disruption of 107 
meteors 


Single-body ablation : 0.3 


The most thorough study of bulk densities through ablation modelling is Kikwaya 
et al (2011). They looked at 107 meteors observed by intensified cameras in 
Ontario, Canada and modelled the deceleration and light curves with an ablation 
model that includes fragmentation (Campbell-Brown & Koschny, 2004). They 
robustly searched the entire parameter space to determine fits to the data. A 
subset of the 92 events not associated with a shower was used. Campbell-Brown et 
al. (2013) inspected 10 meteors using the same model. Orbits and trajectories 
were given for the meteors. With this study, we look at trends in Kg. 


The correlation between K, and bulk density was not as strong as hoped. However, 
a clear relationship between Tisserand Parameter and density was seen. A 
Spearman Rank coefficient of 0.795 was found between density and T,, and 0.441 
between K, and density, indicating a tighter correlation between T, and density. 
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Figure 4: K, against bulk density as found with Figure 5: T, against bulk density as found with 
ablation modelling for 92 meteors. (Kikwaya et al ablation modelling for 92 meteors. (Kikwaya et al | 
PAE) 2011.) 


& 7 a, = Jupiter’s semi-major axis | 
T; = +2 (1 — e*)— cos(i) 
a a] 


a,e,t = semi-major axis, 

eccentricity, inclination of object 
As Jupiter exerts the most dominant influence on comets and asteroids, the 
Jovian Tisserand parameter represents an object’s orbit with respect to the Sun 
and Jupiter. Halley Type Comets have T, below 2, Jupiter Family comets have T, 
between 2 and 3, and most asteroids have T, > 3. However this is not a hard 
cutoff as some comets have T, > 3 and some asteroids have T, < 3. 


When comparing these results to comets and asteroids, one should note that 
these bulk densities are for small particles, ~10° to 10° kg. The macroporosity 
_ of larger objects will change the overall bulk density. Carry (2012) shows 
reasonable porosities between O and 60%. 


P(%) = macroporosity 
p = asteroid bulk density 
Pn —bulk density of associated meteoroid 
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The task of determining meteoroid density distributions is complex and led to the 
investigation described in this report. In reproducing work done by Ceplecha 
(1977, 1988) to correlate K, with bulk density, we examined the modeling of 
high-resolution meteor data from Campbell-Brown et al. (2013) and Kikwaya et 
al. (2011). The correlation between K, and density was not as strong as hoped. 


The Tisserand parameter of a meteoroid is a better indicator of the density than 
the strength proxy, a somewhat surprising result. 


